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ABSTRACT 

A few percent of quasars show strong associated Mg II absorption, with velocities (v ff) lying within a few 
thousand kms" 1 from the quasar systemic redshift. These associated absorption line systems (AALs) are usu- 
ally interpreted as absorbers that are either intrinsic to the quasar and its host, or arising from external galaxies 
clustering around the quasar. Using composite spectra of ~ 1,800 Mgll AAL quasars selected from SDSS 
DR7 at 0.4 < z < 2, we show that quasars with AALs with v s < 1500kms _1 have a prominent excess in 
[On] A3727 emission (detected at > la) at rest relative to the quasar host, compared to unabsorbed quasars. 
We interpret this [O n] excess as due to enhanced star formation in the quasar host. Our results suggest that 
a significant fraction of AALs with v ff < 1500 kms -1 are physically associated with the quasar and its host. 
AAL quasars also have dust reddening lying between normal quasars and the so-called dust-reddened quasars. 
We suggest that the unique properties of AAL quasars can be explained if they are the transitional population 
from heavily dust-reddened quasar to normal quasars in the formation process of quasars and their hosts. This 
scenario predicts a larger fraction of young bulges, disturbed morphologies and interactions of AAL quasar 
hosts compared to normal quasars. The intrinsic link between associated absorbers and quasar hosts opens a 
new window to probe massive galaxy formation and galactic-scale feedback processes, and provides a crucial 
test of the evolutionary picture of quasars. 

Subject headings: black hole physics — galaxies: active — quasars: general — quasars: absorption lines 



1. INTRODUCTION 

In the pop ular Sanders & Mirabel picture for quasar evo- 
lution (e.g.. ISandersetal.l ll988t ISanders & Mirabel 119961 
ICanalizo & Stockton! 1200 U iHopkins et alj 120081) . gas-rich 
major mergers between galaxies trigger intense starbursts in 
ultraluminous infrared galaxies (ULIRGs) and channel am- 
ple fuel into the nuclear region to feed the black hole (e.g., 
lHernquisti ri989). The black hole (BH) growth is initially 
buried in dust until some feedback process (either stellar- 
driven or quasar-driven) blows the dust away and a clas- 
sical quasar emerges. The major merger event also trans- 
forms the host morpholo gy from disks to a spheroid (e.g., 
Toomre & Too m"reiri972l) . The feedback processes invoked 
in this model are also thought to be responsible for shutting 
down star formation in the bulge, and for establishing the tight 
correlations between BH mass and bulge pr operties as ob- 
serve d in local dormant massive galaxies (e.g., Tremaine et al. 
2002). Although such a violent route is not required for trig- 
gering low-to-intermediate luminosity Active Galactic Nu- 
cleus (AGN) activity, this merger-based framework seems 
successful in reproducing an array of observations of lumi- 
nous quasars and massive red ellipticals in the cosmological 
context re.g JWvithe & Loebl2003UHopkins et al.(2 008: Shen 
2009). Nevertheless, many fundamental issues remain unre- 
solved in this model, in particular, the nature and physics of 
these fueling and feedback processes along the evolutionary 
sequence remain largely elusive. 

A useful tool to probe the physical conditions and dynam- 
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ical processes within quasar hosts is intrinsic quasar absorp- 
tion lines. Historically most of the focus has been on broad 
absorption lines (BALs, usually defined as absorption troughs 
broader than 2000 kms -1 ), which are undoubtedly intrinsic to 
the quasar. Here we focus on another class of narrow absorp- 
tion lines called associated absorption lines (AALs), whose 
absorption velocity is close to the systemic velocity of the 
background quasar (z a b ~ Z e m)- AALs are traditionally de- 
fined as narrow absorption troughs (< 500 kms" 1 ) with a ve- 
locity offset v ff within ±3000 kms" 1 of the systemic redshift 
of the quasar 5 . Strong (EW > 0.6 A) low-ionization Mgll 
associated absorbers are pr esent in a few percent of the en- 
tire quasar population (e.g.. lVanden Berk et al. 2008). These 
absorption systems are generally believed to be close to the 
quasar and are explained by either (or a combination) of the 
following scenarios: absor ption by external galaxies cluster- 
ing ar ound the quasar (e.g.. lWevmann et"aLll 1 979t IWild et all 
2008); absorption by halo clouds of the quasar host galaxy 
(e.g., [Heckman et alj 1199 lh : absorptio n by material from a 
starburst wind of the quasar host (e.g., Heck man et al.lll990l) : 
or originating from the vicinity of the black hole (< 200 pc) 
based on partial coverage analysis or variabil ity studies (e.g., 
lHamann et aDl!995b iBarlow & Sargentlll997l) . On the other 
hand, classical intervening absorber systems (with z a b "C Zem) 
are not physically associated with the quasar and are absorp- 
tions due to cosmologically interve ning foreground galaxies 
along the quasa r line-of-sight (LOS) (Bahcall & Spitze rn~969l : 
Bergeron 1986). The strength of these intervening absorbers 
is shown to correlate wit h the associated sta r formation rate, 

measured within 10 kpc dMenard et al.ll201 lh. 

Using composite spectra, IVanden Ber k et al l {2008) 

s Some earlier studies used a larger velocity offset cut ~ 5000 — 
6000km s to define AALs. This is partly due to the fact that the systemic 
redshift for many high-redshift quasars is based on the centroid of broad emis- 
sion lines such as CIV, which has a larger uncertainty compared with the 
systemic redshift determined from Mg II or narrow emission lines. 
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FIG. 1. — An example of our fitting procedure to identify Mgll AALs. 
The black and gray lines show the SDSS spectrum and flux density errors 
around the Mg II region. The red line shows our pseudo-continuum plus 
emission line model fit. The dashed line marks the quasar systemic velocity 
based on the centroid of the broad Mg II line. The cyan line shows the con- 
tinuum+emission line subtracted spectrum and the magenta line shows the 
double-Gaussian fit to the absorption doublet. For this example only N pll iv = 1 
absorption doublet was identified. 

showed that AAL quasars on average have dust extinction 
almost twice that observed in quasars with intervening ab- 
sorbers (e.g., I York et aT]|2006l) . and the associated absorbers 
have a higher ionization state than intervening absorbers. To 
some extent, this supports the idea that AALs are physically 
associated with the quasar (or at least affected by the radiation 
field of the quasar), but a direct link betw een AALs and the 
quasar is still elusive (cf.. lWild et al.ll2008l) . 

Using ~ 1800 quasars with associated Mg II absorbers from 
the SDSS, we will show that there is a link between AALs 
and [On] emission from the quasar host, suggesting that a 
substantial fraction of AALs are intrinsic to the quasar and 
its host. We further argue that the properties of AAL quasars 
are in favor of them being the transitional population in the 
Sanders & Mirabel picture, with signatures of feedback. We 
describe our sample in ^2] and present our main results in §|3] 
In $4] we present an evolutionary picture to interpret our re- 
sults and discuss its implications. 

2. DATA 

We search all the DR7 quasars (ISchneider et al.ll2010l) with 
Mgll coverage (~ 85k quasars) for associated Mgll ab- 
sorption. This restricts our sample to 0.4 < z < 2 (z ~ 
1.2), which extends to lower redshift than the sample in 
IVanden Berk et al] d2008). The unabsorbed continuum plus 
emission line flux i s modeled with a % 2 fit with rejections of 
absorption troughs (IShen et al.ll201 ll) . We refer the reader to 
that paper on the details of our fitting procedure. In short, we 
fit the restframe 2200-3090 A region with a power-law con- 
tinuum plus an iron template, and a set of Gaussians for the 
Mg II line. During the fitting, we mask out 3a pixel outliers 
below the 20 pixel boxcar-smoothed spectrum to minimize 
the effects of narrow absorption troughs. We then subtract 
the pseudo-continuum plus emission line emission from the 
spectrum to get the residual spectrum. Mgn AA2796,2803 ab- 
sorption doublets are then identified on both sides of the Mg II 
emission line by fitting double-Gaussians to the absorption 
doublet. We only keep absorbers that are detected at > 3cr for 
both Gaussian components. We measure restframe equivalent 
widths (EW) of these AALs and uncertainties of EWs are es- 
timated from the Gaussian fits. Fig. [T] shows an example of 
our fitting results. Although this fitting recipe was developed 
to reduce the effects of narrow absorptions on the emission 
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FIG. 2. — Distribution of velocity offset v ff of associated Mg II absorbers. 
Positive (negative) values indicate blueshifted (redshifted) from the systemic 
redshift. The black histogram shows the whole sample while the red one 
represents a subset of quasars with z < 1.4 for which we have emission line 
redshifts based on [O II] or [O III]. Systems with v g > 1500 kms~' are found 
to be similar to classical intervening absorber systems. 
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FIG. 3. — Top: Median composite spectra for different samples, stacked 
in the quasar rest frame and normalized at 3000 A. The three AAL samples 
(intervening-likeiyellow, blueshiftedxyan, and redshifted:magenta) show in- 
creasing reddening relative to the control sample (black), but are substan- 
tially less reddened than the SDSS dust-reddened quasars (with color excess 
A(g—i) > 0.3) shown in red. Bottom: Flux ratios between composite spectra 
of the three AAL samples and the control normal quasars. While the AAL 
samples show (by construction) Mg II absorption with different velocities, an 
excess of [O II] emission is seen for the blueshifted and redshifted samples 
at the quasar systemic redshift. 
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TABLE 1 
The Sample of AALs 



SDSS Designation 


Plate 


Fiber 


MJD 


redshift 


EW 27 96 


EW Err27% 


EW 28 03 


EW Err2803 


v off 


l 'off,SDSS 


v off.HW 


(hhmmss.ss+ddmmss.s) 










(A) 


(A) 


(A) 


(A) 


(kms- 1 ) 


(kms -1 ) 


(kms -1 ) 


000045.77+255106.1 


2822 


339 


54389 


1.4446 


0.61 


0.07 


0.50 


0.07 


1300 


2726 


2069 


000140.70+260425.5 


2822 


322 


54389 


0.7653 


0.48 


0.04 


0.47 


0.04 


47 


54 


113 


000219.64+260029.2 


2822 


367 


54389 


1.2507 


3.07 


0.18 


3.00 


0.14 


-421 


-588 


-289 



NOTE. — The list of AALs used in this paper. In the last three columns we list the absorber velocity offset calculated from the fiducial systemic 
velocity (based on the broad Mg II centroid), the SDSS redshift in the DR7 quasar catalog, and the improved redshift from Hewett & Wild (201 1). 
The full catalog is available in the electronic version. 



line measurements rather than a dedicated recipe for finding 
narrow absorption lines, it does a reasonably good job in iden- 
tifying AALs. With this m ethod we recovered ~ 80% of the 
AALs in the DR3 sample of lVanden Berk et alj d2008l) . where 
most of the "missing" AALs are weak absorptions with other- 
wise normal quasar properties. However, we do not intend to 
quantify the completeness in our AAL selection as functions 
of S/N and absorber strength, as such a task would require 
a more dedicated narrow absorption finder and Monte-Carlo 
simulations, which are unnecessary for the purposes of this 
paper. 

To define an AAL one needs to know the systemic red- 
shift of the quasar. Stellar absorption features are generally 
unavailable due to the overwhelming quasar light. The red- 
shifts based on narrow lines such as [On] and [Olll] are 
generally consiste nt with stellar absorp tion redshifts within 
< 100 kms" 1 (e.g. jRichards et al. 2002). The redshifts based 
on the broad Mg II emission line are consistent with those 
base d on \Q n] or [O ml with a dispersio n of ~ 350kms _1 
(e.g. jRichards et alj2002trShen et alhoiTl their fig. 17). The 
high-ionization broad Civ line is known to be systematically 
blueshifted (~ 700 kms" 1 ) from low-ionization lines such as 
Mgn, with a large dispers i on of ~ 700km s" 1 (e . g., |G askel 
1982; iTvtler & Fan! fl99l iRichards et all 12001 
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201 II) . and hence will bias the systemic redshift determina- 
tion at high redshift. 

Here we adopt the centroid of the broad Mg II emission line 
as the quasar systemic redshift to calculate the absorber ve- 
locity offset, and use 3000km s" 1 as the velocity cut to de- 
fine an AAL. We compared the absorber velocity offsets with 
those calcula ted by adopting t he imp roved redshifts of SDSS 
quasars from iHewett & Wild! (1201 ll) . and found good agree- 
ments with a dispersion of ~ 300kms _1 . This dispersion is 
consistent with the redshift uncertainty ~ 350 kms -1 based 
on the Mg II centroid. Therefore we adopt a nominal uncer- 
tainty of Av ff = 350kms _1 in absorber velocities, which is 
dominated by systemic redshift uncertainties. We finally ar- 
rived at a sample of ~ 1800 quasars with AALs, among which 
~ 10% show mu ltiple AALs in SDSS spectra (consistent with 
IWild et al.l F2008). The sample of AALs used in the following 
analysis is listed in Table Q] Additional information about 
these quasars can be retrieved from the value-ad ded SDSS 
DR7 quasar catalog compiled in lShen et all (1201 ll) . 

Fig.|2]shows the distribution of absorber velocity v g for the 
whole sample in black, and for a subset of z < 1 .4 quasars for 
which we have redshifts from [O n] or [O III] in red. Positive 
values indicate blueshifted relative to the systemic and neg- 
ative values indicate redshifted. We found that the absorber 
velocity distribution using this systemic redshift definition is 
similar (albeit slightly broader due to the less accurate Mg II 
emission line redshifts) to that using the subset of our sam- 



ple with [On] or [O m]-based redshifts. The dispersion in the 
v fi distribution is ~ 600 kms" 1 , larger than the typical uncer- 
tainty Av ff ~ 350km s" 1 arising from systemic redshift un- 
certainties based on Mg II. It is thus dominated by an intrinsic 
dispersion. Our composite spectra (see §|3) also confirm that 
this dispersion is not caused by incorrec t systemic redshifts. 

As earlier studie s alrea dy show (e.g.. We vmann et aill979l 
IVanden Berk et al] [2008), there is a substantial fraction of 
AALs that are redshifted from the systemic velocity. These 
systems could either be truly infalling absorbers, or the emis- 
sion lines used to estimate the systemic redshift have signif- 
icant outflowing velocity. At v ff > 1500 kms" 1 , the distribu- 
tion flattens out and is consistent with the expectation from 
cosmologi cally intervenin g absorbers, as seen in earlier stud- 
ies (e.g.. IWild et al.ll2008l) . We thus argue that these high- 



velocity AALs are mostly classical intervening absorbers, 
which will be confirmed by our composite spectra in §|3] 

We divide our sample into three categories: 1) AAL quasars 
with v D ff > 1500kms _1 (the "intervening-like" sample); 2) 
AAL quasars with < v off < 1500kms _1 (the "blueshifted" 
sample); 3) AAL quasars with -3000 < v ff < Okms" 1 (the 
"redshifted" sample). We have ~ 750 quasars each in the lat- 
ter two samples, and ~ 370 quasars in the former. 

3. COMPOSITE SPECTRA 

To study the average properties of these three AAL quasar 
samples we create high S/N composite spectra in the rest- 
frame of the quasar. SDSS spectra integrate all the light 
received within a 1.5" fiber radius, corresponding to ~ 10 
kpc at z ~ 1, hence they cover a substantial fraction of the 
quasar host. For comparison purposes we construct a con- 
trol sample of more than 9,000 quasars without detected 
AALs, matched in redshift and /-band magnitude to the AAL 
quasars. We create composite spec tra using the method de- 
scribed in Vand en Berk et al.l (120011) . In short, each spectrum 
was shifted to restframe, rebinned onto a common dispersion 
of 1 A per bin, and normalized. The final composite spectrum 
was generated by taking the median (or geometric mean) flux 
density in each bin of the shifted, rebinned, and scaled spec- 
tra. An error spectrum was generated from the 68% semi- 
interquantile range of the flux densities in each bin scaled by 

N S pec, where N spec is the number of spectra contributing to that 
bin. We refer the reader to lVanden Berk et al.l(l2001b for more 
details regarding generating the composite spectra. 

The left panel of Fig. 3 shows median composite spec- 
tra for the above quasar samples. AAL quasars have col- 
ors lying between normal quasars and the so-called "dust- 
reddened" quasars in SDSS (with color excess 6 A(g-i) >0.3; 

6 A(g — i) is a good indicator of in trinsic quasar colors which accounts 
for the band shifting with redshift (e.g., Richards et al. 2003); larger A(g — i) 
values indicate redder colors. 
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FIG. 4. — Left: Continuum-subtracted median composites of the [Oil] emission line, stacked in the quasar rest frame. The black lines show the level of 
emission for the control quasars. We detect a 70 and 53% excess for the blueshifted and redshifted samples respectively. Right: Same quantify shown for the 
[NeV] AA3346.3426 emission line. In this case we do not detect any change in the line flux. 

tion curve, we require E(B-V) ~ 0.12 to achieve the level of 
[O u] emission enhancement relative to the continuum, sub- 
stantially larger than the inferred E(B-V) ^0.03 from the 
composite spectra. 

We now focus on the [On] emission excess. To quantify it 
and estimate its significance, we create continuum-subtracted 
median composites around the [On] emission line. This is 
done by subtracting a running median filter with a width of 
20 pixels and by masking the region 3710 < A < 3750 A to 
preserve the emission line. Within this region the continuum 
estimate is interpolated. The results are shown in Fig. |4] The 
centroid of the excess [On] emission is within ~ 150kms _1 
from the systemic velocity (Fig. [4j. Considering the typical 
redshift uncertainties of individual objects used in the stack- 
ing this is consistent with the systemic velocity. We measure 
the line flux by integrating the flux density over the emis- 
sion line profile. The line flux error is estimated by boot- 
strapping the sample of reference quasars and repeating the 
procedure 100 times. This allows us to include sample vari- 
ance in the error estimation. The [O n] excess is detected at 
8rj (la) for the blueshifted (redshifted) sample, correspond- 
ing to an enhancement of 70% (53%) of the line flux of the 
control quasars. The enhancement increases to 94% (63%) if 
we use mean composites. For comparison, we measure the 
properties of the [Ne v] AA3346,3426, which is entirely ex- 
cited by AGN given its high ionization potential. We follow 
the same procedure and the results are presented in the right 
panel of Fig. [4] In this case, we do not observe any signif- 
icant excess in the [Nev] AA3346,3426. Therefore the AGN 
ionizing source of AAL quasars is similar to that of normal 
quasars, and we conclude that the excess [On] emission is 
due to different host properties. Interestingly, no [O u] excess 
is seen in the composite spectrum of ~ 300 q uasars with broa d 
Mg II absorption lines selected from SDSS (IShen et al1l2011l) . 
which is expected if broad absorption lines (BALs) are ubiq- 
uitous among quasars and the appearance of BALQSOs only 
depends on LOS. This difference reinforces our conclusion 
that the hosts of AAL quasars are intrinsically different. 



Richards et al. 2003). The reddening of AAL quasars rel- 
ative to the control quasars is well described by a SMC- 
like extinction curve, wit h E(B-V) ~ 0.03, consistent with 
IVanden Berk et al] (|2008i) . The right panel shows flux ratios 
of the AAL quasar composites to that of the control quasars. 
The "redshifted" and "blueshifted" samples show a prominent 
excess of narrow [O n] emission, which is absent in the case of 
the "intervening-like" sample. The lack of [O n] emission ex- 
cess, the similar amount of reddening to classical intervening 
absorption systems (e.g., York et al. 2006), and the fact that 
they join the plateau of the velocity distribution (e.g., Fig.|2j, 
suggest that AALs with v g > 1500kms _1 are in fact mostly 
classical intervening absorbers 7 . We note that the large-scale 
fluctuations seen in the composite ratios are due to correlated 
noise arising from variance in the shape of quasar continua. 

The correlation between the presence of absorbers with 
v D ff < 1500 kms" 1 and [On] emission is of great interest. 
While the composite Mg II absorption lines are, by construc- 
tion, offset by hundreds of kms" 1 with respect to the quasars, 
the stacked excess [O n] emission is found more or less at the 
quasar systemic velocity (see below). Thus the excess [On] 
emission originates from material lying at the quasar systemic 
velocity. The central engine, i.e., the black hole radiation, ap- 
pears to be similar in AAL quasars and in normal quasars: 
other than the reddening and Mg II absorption, no difference 
is seen in the continuum and broad emission lines between the 
two quasar populations. This suggests that dust-reddening is 
the explanation for the color difference seen in AAL quasars 
and normal quasars. One might ask if this dust reddening 
might also cause the apparent [On] excess seen in the flux 
ratio plot in Fig. [3] Dust reddening could occur on spatial 
scales much smaller than the [O n] emission region but much 
larger than the continuum plus broad line region. In this case 
the continuum is attenuated while the [O n] emission is not, 
causing an apparent enhancement of [On] strength relative 
to the underlying continuum. Assuming an SMC-like extinc- 

7 This statement is for the low-z Mg II AALs, and may not hold true for 
the high-z and high-ionization AALs (such as CIV). 
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4. DISCUSSION 

The distinct host properties of AAL quasars in terms of ex- 
cess [On] emission (with absorber v fr < 1500 kms" 1 ) suggest 
that a significant fraction of these associated absorbers must 
be physically associated with the quasar and its host, and are 
not from disconnected external galaxies 8 . Because the prop- 
erties of the active SMBH are similar for AAL quasars and for 
normal quasars, we expect that the amount of [O n] emission 
contributed by the quasar itself should be similar as well. We 
thus interpret the excess [O n] emission as due to enhanced 
star formation within the AAL quasar host galaxy, as [O u] 
is a good indicator of star formation rate in regions with low 
ionizatio n para meters, and has been applied for quasar hosts 
(e.g.. lHoll2005l) after subtractin g the [O n] contribu tion from 
the quasar itself. Interestingly, iMenard et al.l (1201 ll) showed 
a similar connection between Mg II absorbers and star for- 
mation in intervening galaxy absorbers. The measured ex- 
cess of [O n] luminosity of the AAL quasar hosts is found to 
be L = 1.9 and 1.2 x 10 4I ergs _1 for th e blueshifted and red - 
shifted samples respectively. Following Me nard et al.l (i20 11). 
we convert these values to star formation rates using the aver- 
ag e scaling coefficie nt for galaxies with 0.75 < z < 1 .45 given 
bv lZhu et alJ J2009b . We then find AAL quasar hosts to have 
a median excess star formation rate of ^7 (5) M Q yr" 1 for the 
blueshifted (redshifted) sample. 

A less likely interpretation for the [O n] excess would in- 
voke different ISM conditions or element abundance for AAL 
quasar hosts. Either way, our results show that AAL quasars 
represent a distinct population and may indicate a special evo- 
lutionary phase. The clear absence of [On] emission at the 
redshift of the associated absorbers in the composite spec- 
tra disfavors major mergers (in which case enhanced star for- 
mation for the companion would also have been seen). Our 
results thus imply that associated Mg II absorbers could ei- 
ther be minor mergers, or galactic to halo-scale gas outflows 
and infall associated with the quasar host galaxy. We still 
expect some AALs from the clustering of galaxies around 
quasars (e.g.. IWild et al.l i2008). Quantifying their contribu- 
tion requires better knowledge on the strength of the galaxy - 
quasar cross-correlation function and pairwise velocity dis- 
persion (for the specific type of absor bers), as well as the 
size o f the quasar proximity zone (e.g jHennawi & Prochaskal 
2007), and is beyond the scope of this paper. 

A number of studies have revealed that dust-reddened 
quasars have different properties from classical quasars: heav- 
ily reddened (with E(B-V) > 0.5) type 1 quasars, selected 
in the near infrared and radio, usually show strong [On] 
emission (e.g., Glikman et al. 2007). In addition, these 
systems display a large fraction (approaching unity) of dis- 
turbed morphologies and interaction (e.g., Urrutia et al. 2008). 
Dust reddened quasars from the SDSS (with E(B-V) ~ 0.1) 
also show an excess of [On] emission compared to classi- 
cal quasars (Richards et al. 2003). This is in line with the 
Sanders & Mirabel picture that dust-reddened quasars are 
in the early stage of evolution after the merger and have 
higher star formation rate and more dust obscuration in the 
nuclei. It is interesting to note that the reddest quasars also 
show a higher oc currence of associated absorbers than nor- 
mal quasars (e.g., lAldcroft et al.lll994t iRichards etaD 120031: 
HopkinsetaDll004), but the bulk of AAL quasars have much 
bluer colors (Fig. 3). 

8 Although these external galaxies may still be affected by the ionizing 
flux from the quasar. 



AAL quasars thus have properties intermediate between 
dust-reddened and classical quasars, in terms of both [O u] 
emission and dust reddening. These findings are consistent 
with AAL quasars being more evolved than the dust-reddened 
quasars following a gas-rich merger. The obscuring material 
in the immediate circumnuclear region has largely been dis- 
persed, but the host star formation rate is still higher than nor- 
mal quasars. Alternatively, AAL quasars could be similar to 
dust-reddened quasars, but are seen along a less obscured line 
of sight. In either scenario, we expect AAL quasars to show a 
higher fraction of less-developed bulges, disturbed morpholo- 
gies and interactions than classical quasars. 

An archival search in the HST database resulted 
in three AAL quasars in our sample that have deep 
(> lOmin) imaging data at z ~ 0.6, all observed in 
serendipitous programs and the data is publicly avail- 
able. Two of them, SDSSJ085215.65+492040.8 and 
SDSSJ130952.89+01 1950.6, were observed with WFC3-IR 
(~ 17 min and ~ 13 min) and WFC3-UVIS (~ 28 min and 
- 25 min). The third object, SDSSJ093210.96+433813.1, 
was observed with WFC3-IR with a total exposure time of 
~ 40 min. In Fig. [5] we show the retrieved WFC3-IR im- 
ages for these three objects, and their SDSS spectra around 
the Mgll region. J0932 and J 1309 show clear tidal features 
and disturbed host morphologies even before subtraction of 
the quasar light, while J0852 appears to be in a small in- 
teracting group where the immediate neighbor to the lower- 
right has a SDSS photometric redshift consistent with the 
quasar redshift. Sluse et al. (2007) reported a serendipi- 
tously discovered AAL system in the gravitationally lensed 
quasar RXS Jl 131-1231, which also shows evidence of inter- 
actions and disturbed morpho logy in the reconstructed host 
image (Claeske ns et alJ 12006). Thus the frequency of dis- 
turbed morphologies and close compani ons of AAL quasars 
is much higher than e volved quasars (e.g.. | Bahcall et al.ll!997t. 
Dunlo pet al.l 120031 : IVeilleux et alJ l2009h . and is more in 
line with the result s for heavily dust-reddened quasars (e.g., 
lUrrutia et al.ll2008l) . This is fully consistent with our evolu- 
tionary scenario, but more data are needed to draw firm con- 
clusions. We defer a detailed analysis of the HST data and 
follow-up studies of these systems to future work. 

Even though our results suggest an intrinsic origin for the 
bulk of the AALs, this is a statistical statement and there are 
still cases where the AALs are of external origins. Moreover, 
the nature of these intrinsic AALs is still unknown. Never- 
theless, they provide a unique window to probe the physical 
conditions and dynamical processes during the evolution of 
quasars and their massive hosts, and deserve in-depth studies, 
either in a statistical manner or for individual systems. More 
than half of these associated absorbers show blueshifted ve- 
locity, among which some could be feedback (either stellar- 
or quasar-driven) at work, and possibly leading to a suppres- 
sion of the host star formation. The typical velocity of these 
absorbers (i.e., a few hundred kms" 1 ) is much smaller than 
those seen in BALs and intrinsic high velocity narrow absorp- 
tion lines (NALs), hence these AAL absorbers are likely to 
be on scales much larger than nuclear scales 9 , and thus are 
more likely to impact the host galaxy. Gravitational inflows 
may also be present during the early stage of quasar and host 

9 Associated absorptions are therefore expected to show little variability on 
multi-year timescales compared with BALs and high- velocity intrinsic NALs, 
which can be tested with multi-epoch spectroscopy of these absoiption line 
systems (Shen et al., in preparation). 
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FIG. 5. — Three AAL quasars in our sample with deep HST images. Upper: WFC3-IR images with asinh scaling. From left to right: J0852+4920 (F125W), 
J0932+4338 (Fl 10W) and J1309+01 19 (F125W). The quasars are located near the centers of the images. Bottom: The Mg II region of the optical spectrum from 
SDSS for the three objects. 
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formation, either in the form of accreted gas or as part of an 
infalling galaxy companion. While intrinsic absorbers with 
infalling velocities greater than several hundred kms" 1 are 
likely to be in the vicinity of the BH, the nature of the more 
slowly moving absorbers needs to be explored. In particular, 
infalli ng absorbers are e xpected in galaxy formation models 
(e.g., iKeres etal J 120051: iDekel et alj I l2009h but almost never 
seen in galaxy spectra (Steidel et al. 2010). The possibil- 
ity that some of these infalling absorbers might be associated 
with cold flows will be explored in future work. 
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